Spin hydrodynamics in the S = 1/2 anisotropic Heisenberg chain 
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We study the finite-temperature dynamical spin susceptibility of the one-dimensional (generalized) 
anisotropic Heisenberg model within the hydrodynamic regime of small wave vectors and frequencies. Numer- 
ical results are analyzed using the memory function formalism with the central quantity being the spin-current 
decay rate y(q, oj). It is shown that in a generic nonintegrable model the decay rate is finite in the hydrodynamic 
limit, consistent with normal spin diffusion modes. On the other hand, in the gapless integrable model within 
the XY regime of anisotropy A < 1 the behavior is anomalous with vanishing 7(5, ui — 0) oc \q\, in agreement 
with dissipationless uniform transport. Furthermore, in the integrable system the finite-temperature q = dy- 
namical conductivity a(q = 0, uj) reveals besides the dissipationless component a regular part with vanishing 

0"reg(<? = 0, U) — > 0) — > 0. 

PACS numbers: 05.60.Gg, 71.27.+a, 75.10.Pq 



I. INTRODUCTION 

The S = 1/2 anisotropic Heisenberg model (AHM) on a 
one-dimensional (ID) chain is one of the prominent models 
describing the physics of strongly interacting fermions on a 
lattice. The model is also well realized in several novel mate- 
rials, e.g., in the ID Mott-Hubbard insulator the relevant low- 
energy degrees of freedom remain spin excitations, and the 
closest realizations correspond to the isotropic case A = 1. 
The integrability of the model via the Bethe Ansatz has so far 
led to a number of exact results regarding the spin dynamics 
at T = 0, 1 as well as some thermodynamic results at T > 0. 2 
On the other hand, it has been also realized that just the inte- 
grability itself, of this particular model or more generally, can 
be the origin of the anomalous behavior of transport quantities 
and low-energy dynamics. 3-6 

A specific consequence of model integrability on a chain 
of L sites is the existence of a macroscopic number of con- 
served local quantities and related operators Q n , n = 1 , . . . , L 
commuting with the Hamiltonian [Q n , H] = and with each 
other [Qn,Qm] = 0. For the ID AHM a nontrivial exam- 
ple is Q3 = J E representing the energy current and leading 
directly to its non-decaying behavior 7 and within the linear 
response to infinite thermal conductivity k(T) at any T > 0, 
being also the explanation for the very large spin contribution 
to heat conductivity in spin-chain materials. 8 

We concentrate in this paper rather on spin dynamics and 
response, as evidenced by the dynamical spin susceptibil- 
ity x(g,w), with the emphasis on the hydrodynamic regime 
of small (q, oj) — > 0. The behavior is also expected to be 
anomalous since it has been shown that spin conductivity 
a(q = 0,u>) has a ballistic (dissipationless) contribution - 
spin stiffness D(T) > - again at any T > within the 
gapless XY regime A < l. 6 Although the corresponding spin 
current Jq is not a conserved quantity, the ballistic component 
has been well established via several connections: (a) the spin 
current is closely related to many-body (MB) level dynamics 
€n{4") w i tn respect to flux (f> 9 and their independent-like char- 
acter in an integrable MB model 3 ' 10 in contrast to level re- 
pulsion in a generic nonintegrable system, (b) via the limited 



decay of correlation functions due the overlap with conserved 
quantities 7 , ( J§ Q3) ^ 0, at least for a partly polarized sys- 
tem with nonzero magnetization, s 7^ 0, (c) the overlap with a 
steady state solution 1 1 which - unlike other local conservation 
laws - is not orthogonal to the spin current in the s = sector. 
While the behavior is evidently different 12 but as well highly 
nontrivial in the gapped Ising-like regime s = 0, A > 1 
where results seem to favor D(T) = fj, 10 ' 11 ' 13 there seems 
to be most controversy on the strict vicinity of the isotropic 
point A ~ 1, s ~ o. 6 ' 10,14 There is much less known about the 
T > spin dynamics extended to finite wave vectors q > as 
well as frequencies oj > 0, which has been considered so far 
in finite systems using exact diagonalization 1516 and quantum 
Monte Carlo 1718 (QMC). While for a generic nonintegrable 
spin system the hydrodynamic regime should reveal diffusive 
behavior for (q, u) — > 0, in the gapless integrable AHM a sin- 
gular approach at q — > is expected. In the latter case one 
of the open questions is the possible coexistence of a ballistic 
and diffusive processes at q — s- 0. 19 ' 20 Clearly, the central chal- 
lenge is the design of a proper phenomenological description 
of the T > hydrodynamic response of the gapless integrable 
as well as weakly perturbed nonintegrable spin system. 

In this paper we analyze spin hydrodynamics in the inte- 
grable and nonintegrable AHM by the numerical calculation 
of T > dynamical spin correlations on finite chains. Nu- 
merical data are used as the input to the memory-function rep- 
resentation where the crucial result is the spin-current decay 
rate 7(9, w) at T > 0. The latter quantity should be constant 
7(g, oj) = 70 within the diffusion regime, as is shown to be 
for the nonintegrable case. For the integrable gapless AHM 
within the XY regime, however, our results indicate vanishing 
j(q -> 0,oj = 0) = as well as 7(17 = 0,w -> 0) = 0. This 
is well consistent with the dissipationless transport at q = 
but puts also restrictions on possible proposed theoretical sce- 
narios. 

The paper is organized as follows: In Sec. II we present the 
model, spin correlations and the memory-function formalism 
for the dynamical spin susceptibility. Section III is devoted 
to the presentation of numerical methods. We focus also on 
the evaluation of the decay rate j(q,uj) from dynamical spin 
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correlation functions. Results are presented in Sec. IV and 
Sec. V. First we investigate transport properties in a nonin- 
tegrable case with a next-nearest neighbor interaction. Next 
we repeat our investigation for high- and low-T behavior of 
q dependent j(q,uj) in the case of the integrable AHM. We 
also extend our analysis for the case of finite magnetization 
s 7^ 0. Finally, we focus in Sec. VI on the low-w behavior of 
the uniform dynamical conductivity. 



Applying Mori-Zwanzig reduction 21,22 we get 



$(?,«) = 



(6) 



UJ + tq 2 a(q,uj)/Xq ' 
Here, <j(q, uj) is the projected dynamical spin conductivity, 

a(q,uj) - »(j;|(cj - QCQ)~ l \r Q ) , (7) 



II. MEMORY-FUNCTION ANALYSIS 

In this paper we study the anisotropic Heisenberg (XXZ) 
model on a chain with L sites and periodic boundary condi- 
tions (PBC), 



where the spin current is defined by the conservation law S q — 
iqJq (written at small momentum q <C 1 ) and the orthogonal 
projection reads Q = 1 — \S*)(S% \/x q - The spin diffusion 
constant (if it exists) is equal to T> = lim g ^o &(q, w = 0)/x°- 
On the second level we can represent 



H = JJ2(S*S\ 



+i 



J i °i+l 



AS?S, 



IS.2SIS, 



(1) 

where Sf (a = x,y,z) are spin S — 1/2 operators at site 
i and A represents the anisotropy. We allow also for a next- 
nearest neighbor zz-interaction with A2 7^ breaking the in- 
tegrability of the model. It should be reminded that the Hamil- 
tonian (1) can be mapped on a t-V-W model of interacting 
spinless fermions with hopping t = J/2 and inter-site inter- 
actions V = J A, W = JA2. In the fermionic representation, 
spin transport corresponds to charge transport. We further on 
use h = ks = 1 as well as J = 1 as the unit of energy. 

Our aim is to investigate the dynamical spin susceptibility 
X(q,u) given by 



(2) 



X(5,«)=* / dte^{[S z _ q {tlSl}) 



with 



si 



(3) 



where due to PBC q = 2nk/L and (. . .) denotes the thermo- 
dynamic average at temperature T. Note that the dynamical 
spin structure factor is related to the imaginary part of Eq. (2) 
by S(q, uj) = tt[1 - exp(-u>/T)} X "(q, u). 

While the results for x(Qi^) are obtained numerically 
for finite systems, we use for further analysis the memory- 
function (MF) formalism 2 1-23 with the central quantity being 
the relaxation function 



<£>(q, uj) 



X(Q,v) - Xg 



= (S*\(C-co)- 1 \S*), (4) 



where \q = x(<l> 0) = (55 \Sg) is the static spin susceptibil- 
ity, C is the Liouville operator, LA = [H, A], and within the 
MF formalism the scalar product between operators at T > 
and corresponding = 1 /T is defined by 



B) = ^J At{A^B{it)) . (5) 



a(q,uj) 



lj + M(q, uS) 



(8) 



where x° q = {S z q \S z q )/q 2 = {J q \J q ), such that 



-v° 



(9) 



u: + M (q, uj) 



The memory function M (q, uj), and in particular its imagi- 
nary part j(q,u)) = M"(q,uj) representing the spin-current 
decay rate, is the central quantity of our analysis. In the 
generic case of normal transport it should be well behaved 
in the hydrodynamical regime (q, uj) — > with a finite value 
70 = j(q — > 0,uj — > 0). This is clearly not the case 
for the dissipationless (integrable) case where one can expect 
j(q -> 0,w = 0) -> 0. 

On the other hand, Xq an d are t> om finite in the limit 
q 0. Moreover, both are even q independent at high T ^> 
(j3 — > 0), as can be easily shown via the high-T expansion, 



Xq 



= /3 (1/4 - s 2 ) + 0{(i 2 
/3^(l/4-,s 2 )+0(/3 2 



(10) 



where s = (Nf — N^)/2L is the magnetization and (A^) 
is number of up (down) spins. 

It should be pointed out that, due to the projection Q in 
Eq. (7), a(q, u>) is not equal to the "standard" spin conductiv- 
ity cr(q,uj), 



<j(q,Uj) = l(J q \{uJ - C) 1 \J q ) . 



(11) 



Note that the latter quantity is directly related to the relax- 
ation function a'(q,uj) = uj 2 Q"(q,uj)/q 2 and thus differs 
essentially from a'(q,uj) for q > 0, in particular in the 
regime uj — > 0. On the other hand, for strictly q = (in 
a finite system with PBC) both quantities should be equal, 
a(q = 0, uj) = a(q = 0, uj), but in this case one has to calcu- 
late numerically er(0, uj) and corresponding M(q = 0, uj) sep- 
arately [not from <3?"(0,u;)] since S q=0 is a conserved quan- 
tity. 
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III. NUMERICAL METHODS 

In this paper we present results obtained using different nu- 
merical methods: 

(a) For high T > 1 studies we use the Microcanonical Lanc- 
zos Method 24 26 (MCLM). The choice is motivated by the fact 
that in the thermodynamic limit the microcanonical ensemble 
should yield the same results as the canonical one. For a large 
enough system L and temperature T, dynamical autocorrela- 
tions can be evaluated with respect to a single wave function 
characterized by the energy uncertainty 



<5e=V«*l(#-A) 2 |*» : 



(12) 



where the parameter A = (H) determines the temperature for 
which | "J) is a relevant representative. Such \*f?) can be gen- 
erated via a first Lanczos procedure using (H — A) 2 instead 
of H. The dynamical correlation functions are then calcu- 
lated using the standard Lanczos method, where the modified 
wave function $ = S*\*S!) is the starting point for the second 
Lanczos iteration. Reachable finite-size systems, L = 28 for 
magnetization s = and L = 36 for s = 1/4, have very high 
density of states for high T, hence statistical fluctuations are 
effectively smoothed out. This is in contrast to low-T prop- 
erties, dominated by a small number of low-lying MB states, 
therefore a different numerical method should be applied, 
(b) For low T < 1 we choose the Finite-Temperature Lanczos 
Method 25 - 26 (FTLM) which is reliable for T > T fs where the 
finite size temperature is typically T ts ~ 0.5 in the AHM at 
available chain lengths L. To reduce statistical fluctuations at 
low T, one can introduce additional sampling over initial ran- 
dom vectors and increase the number of Lanczos steps. Due to 
the memory requirement and the CPU time, reachable system 
sizes are L = 26 for s = and L = 32 for s = 1/4. 

Within the framework of linear response theory we numer- 
ically calculate the imaginary part of the dynamical suscepti- 
bility x"(<Z,k-0 and for q = the spin dynamical conductiv- 
ity cr'(0, w). The MF 7(5,0;) can be calculated directly from 
Eq. (9), 



7(5,0;) 



(13) 



However, the denominator of this equation is highly influ- 
enced by the Kramers-Kronig transformation, e.g., by a small 
imaginary process e added in the numerical realization of the 
Kramers-Kronig relation, 



$(q,u>) 



.1 



uj — ui' + it 



(14) 



To resolve this problem it is better to perform the Hilbert 
transform on x"(<7, w) or cr'(0, w) and use one of the following 
expressions for the MF, 

^a'(q = 0,cu) 



7(9 = 0,w) = 
7(5,0;) 



K<7 = 0,oj)| 2 ' 

_ q 2 Xqx"(<i,u) 
u\x(q,u)\ 2 



(15) 
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Figure 1. (Color online) 7(5,0; = 0) vs. q (in real units) for A = 
0.5, s — 0, as calculated for a integrable system of L = 22, 24, 26, 
and 28 sites. 



which comes from simple evaluation of Eq. (8) and Eq. (9) 
respectively. Finally, we can calculate >P q using one of the 
sum rules 



do; x"(g,w) 
— w 2 



-o° = / ^V( 9 = 0, W ), 



(16) 

where the integration goes over the whole, numerically ob- 
tained MB spectrum, i.e., oj spatl ~ L J. It is also worth men- 
tioning that Xq an d Xq are weakly q dependent in the hydro- 
dynamic limit, as evident for (3 — > 0, cf. Eq. (10). 

As a final remark of this section we show in Fig. 1 result 
for the decay rate 7(5, uj — 0), evaluated for different system 
sizes L = 22 - 28. We choose here A = 0.5 and A 2 = 
case, but this behavior is generic for "conducting" regime. It 
is evident that oj = limit of decay rate revel only system- 
atic size (q) dependence, without sizeable statistical error or 
deviations. 



IV. NONINTEGRABLE CASE 

First, let us address the question of spin transport in the 
generic "normal" case, i.e., nonintegrable case as introduced 
by A2 7^ 0. In Fig. 2 we present characteristic numerical re- 
sults within the XY regime A = 0.5 at high T, as obtained via 
the MCLM method for a chain of L = 28 sites and different 
allowed smallest q„ = n(2ir/L), n = 0, . . . , 3. 

Shown are the dynamical relaxation function $"(q,o;) = 
x" (q, oj) /uj, and the extracted projected dynamical conductiv- 
ity a'(q, uj), Eq. (8), as well as the decay rate 7(5, oj), Eq. (15). 
Results obtained here for a substantial integrability breaking 
term A2 = 0.6 can be easily interpreted with normal diffu- 
sion behavior. We note that the projected conductivity a(q, 10) 
reveals a nearly q independent Lorentzian form, which is then 
reflected in 7(5,0;), being effectively constant in both uj and 
q. In particular, it is evident that the q = result 7(5 = 0, uj) 
obtained directly from a(q = 0,o;) is essentially the same as 
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Figure 2. (Color online) High-T (a) spin relaxation function 
$"(<?, w), (b) projected conductivity a(q,uj) and (c) spin-current 
decay rate y(q, oj) for a nonintegrable AHM at s — 0, T — 10, 
A = 0.5, A2 = 0.6, as calculated for a system of L — 28 sites. 
With dots we present the comparison with a q = perturbation the- 
ory. 
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Figure 3. (Color online) Decay rate 7(3, oj = 0) vs. momentum 
q for the case as in Fig. 2 with A2 = 0.6, but for different A = 
0.5, 0.8, 1.0, evaluated for L = 28. Horizontal lines represent q — 
results from the perturbation theory. 



7(51,0;) extracted from <&"(qi,iS). Deviations for q > q\ are 
understandable since for our restricted chain lengths the latter 
g's are actually not very small. 

As a final numerical result of this section we present in 
Fig. 3 the finite g-scaling of the d.c. rate j(q,oj — 0) for 



three different values of the anisotropy A = 0.5, 0.8, 1 in 
the case of non-integrability A 2 = 0.6 of the model (1). As 
already evident from Fig. 2, the variation with q is modest and 
7(q,oj = 0) is finite in the limit q — > 0. Furthermore, the 
value agrees very well with a perturbation theory at q = on 
the basis of Ref. 27, on which we comment in more detail for 
the remainder of this section. 

To this end, it is convenient to turn to the time domain for 
the moment and to consider an integro-differential equation of 
the form 



d t a'(q = 0,t) 



J drK(t 



T )a'(q = 0,T) (17) 



describing the decay of the dynamical conductivity a'(q = 
0, t) in time and involving a time-convolution with a memory 
kernel K{t). Such an equation, and particularly the memory 
kernel, can be obtained by an application of the Nakajima- 
Zwanzig (NZ) projection operator technique. 28 For the appli- 
cation of NZ an unperturbed Hamiltonian Hq has to be cho- 
sen, where a natural choice is given by the XX model, i.e., 
H = H(A = A 2 = 0). Then the interaction V = H - H 
may play the role of a small perturbation and NZ allows to ob- 
tain the memory kernel K (t) as a lowest-order truncation of a 
systematic series expansion (in even powers of the perturba- 
tion strength involving A and A2). The somewhat lengthy 
and subtle calculation of the lowest-order truncation K2(t) 
has already been undertaken and is given in Eq. (30) of Ref. 27 
by K2(t) = dtR^it). While the final expression for K2(t) 
still needs to be evaluated numerically, the evaluation can be 
done for several thousands sites, e.g., L = 2000 as chosen 
here. Thus, numerically integrating the above Eq. (17) by 
the use of K(t) w K 2 {t) for L = 2000 eventually leads 
to a lowest-order prediction in the thermodynamic limit for 
cr'(q = 0, t), and after Fourier transforming, for a'(q = 0, u>) 
as well. The obvious agreement in Fig. 2 also includes small 
deviations from a strict Lorentzian, which is a non-Markovian 
effect since the memory kernel K2 it) cannot be considered 
as a mere 5 function on the characteristic scale set by the 
decay time of a'(q = 0,t). Notably, K 2 {t) directly yields 
a lowest-order prediction on the frequency-dependent decay 
rate 7(9 = 0, oj) as well, via the relation 



j(q = 0,oj)fa / dt cos(ujt)K 2 {t) , 



(18) 



with an as convincing agreement in Figs. 2 and 3. It is fur- 
ther worth mentioning that the success of NZ relies not only 
on the limit of small perturbations but also crucially on the 
non-integrability, allowing to assume vanishing higher-order 
contributions (equivalent to neglecting the coupling to other 
observables in the NZ equation). In the case of integrability, or 
in the limit of strong perturbations, the incorporation of higher 
order contributions is indispensable and also other variants of 
projection operator techniques may become convenient 27 , not 



containing a memory kernel at all. 




Figure 5. (Color online) Low-T (a) &"(q, uj) and (b) 7(5, uj) for the 
integrable case as calculated for L = 26, s = 0, A = 0.5, and 
T = 0.7. Points represent the analytical A = results for the same 
parameters and q. 




Figure 4. (Color online) High-T (a) $"(q,u), (b) a(q,uj) and (c) 
y(q, uj) for the integrable case at s = 0, A = 0.5, evaluated for 
L — 28 and T = 10. Points represent the analytical A = results 
for the same parameters and q. 



V. INTEGRABLE MODEL 

In the following we analyze the spin dynamics within the 
integrable AHM (XXZ) model as realized for A 2 = 0. We 
concentrate here on the gapless ("conducting") regime, i.e. on 
the XY regime A < 1 at zero magnetization s = 0, but ex- 
tending also to A < 1 for s ^ 0. Note that a different behavior 
can emerge in the "insulating" (spin-gap) phase at A > 1 and 
s = (even at T > 0). 



A. XX model 

In the particular case A = the AHM (i.e. the XX model) 
can be mapped on the tight-binding ID model of noninter- 
acting spinless fermions. In this case the calculation of the 
dynamical susceptibility, Eq. (2), reduced to the well known 
Lindhard formula and for q — > the relaxation function 
$"(<j, ui) can be evaluated analytically for arbitrary T > 0, 







Using the above expression and relation (13) one can express 
directly the decay rate j(q, uj), and in particular 



1 XX (q,uj = 0) = 



[2 X »cosh(/3/2)] 



(20) 



4cosri {K/2q)t 



(19) 



Since the prefactor at any T > is a constant at q — » 0, we 
note a characteristic linear variation j xx oc \q\ and its (non- 
analytic) vanishing for q — > 0. As we see in the following, the 
XX model can serve as the guideline and explanation for the 
observed behavior for more general A > 0. 



B. XXZ model 

Next let us focus on the case of a more general AHM with 
A > (XXZ model) which corresponds, via the Jordan- 
Wigner transformation, to a model of interacting spinless 
fermions and analytical results are not available for T > 
and q > 0. We are therefore restricted to the numerical cal- 
culation of $(q, u>). As described in Sec. Ill, two approaches 
are used: (a) the MCLM for high enough T > 1 and (b) the 
FTLM for lower T. 

In Fig. 4 we first present the high-T results for A = 0.5 
and s = 0. Again we show first $>"(q, u), then the extracted 
a'(q,uj), and finally 7(5,0;). The difference to a noninte- 
grable case is quite evident. Both a'(q,ui) as well as j(q,ui) 
are strongly q dependent, in particular at low u>, most pro- 
nounced for q = 0. Here, tj'(q = 0, uj) = <r'(q = 0,uj) ex- 
hibits a dissipationless component 2ttD5(uj) reflected in the 
vanishing j(q = 0, uj — > 0) — > 0. 

To stress the correspondence of the results to the XX case 
we present in Figs. 4 and 5 results for uj) and j(q, uj) 
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Figure 6. (Color online) 7(9, ui = 0) vs. q for s = and different 
A < 1, (a) low T = 0.7 and L = 26, (b) T » 1 and L = 28. 
Gray-dashed lines are guides to the eye. 
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Figure 7. (Color online) "f(q,ui = 0) vs. <jr for s = 1/4 and different 
A < 1, (a) low T = 0.7 and L = 32, (b) T » 1 and L = 36. 
Gray-dashed lines are guides to the eye. 



in a direct comparison with the analytical A = expressions 
Eq. (19) with the same parameters, same q, as well as T = 10 
and T = 0.7, respectively. Similarities and differences to 
the XX case are quite visible: (a) singularities and cutoffs 
in $"(q,w) at ui c (q) = cq persist at A > as maxima di- 
minishing with increasing A, (b) the qualitative behavior of 
7(5, ui) is quite similar for the A = and A = 0.5 cases for 
lower u) < w c , (c) for A = decay rate vanishes above the 
threshold, i.e., y[q,ui > cj c (q)] = 0, for A > on the other 
hand j[q, co > oj c (q)] becomes q independent but nonzero in 
a broad u> range. 

Next we focus on the finite ^-scaling of j(q,uj). In Fig. 6 
we present -f(q,u) = 0) for the zero-magnetization s = 
as a function of q for various values of the anisotropy (A = 
0, 0.3, 0.5, 0.8, 1) and for two temperatures T = 0.7 and 
T = 10. As we mention before, the q = and g ^ results 
are calculated from different (numerically obtained) quanti- 
ties, dynamical conductivity <r'(q = 0, w) and spin suscepti- 
bility x"(Qi w )' respectively. It is indicative that j(q, 10 = 0) 
within the XY regime is linear function of q. Also our results 
show vanishing 7(5 — > 0, ui = 0) — > for A < 1. This type 
of behavior is well consistent with dissipationless (ballistic) 
transport for zero wave vector q = 0. It is, however, also vis- 
ible that the isotropic case A = 1 does not follow the simple 
7(5, lj = 0) oc q behavior (similar deviations can be observed 
also for A = 0.8 and low T = 0.7 which could be also due 
to finite-size or FTLM numerical uncertainty). As discussed 
later the conclusion could be that -f(q — >• 0, u> — 0) = 7$ re- 
mains finite 29 or even more that the scaling 7(5, lj = 0) oc q a 
is different 30 . 

We also repeat the same analysis for the case of finite 
magnetization s ^ . Figure 7 depicts the q-scaling for 
s = 1/4. Since we are away from the singular isotropic point 
s = 0, A = 1 and D(T) > for all A, 7 here the A = 1 case 
does not show any deviations from a general rule. 

Finally let us focus on the temperature dependence of 
•y(q, lj). As mentioned, the FTLM is reliable for T > Tf s ~ 




Figure 8. (Color online) ^/(q,oj = 0) vs. T for A = 0.5, s — 
0, as calculated for a system of L — 26 sites. Inset: For A = 1 
the comparison between 7(51,0; = 0) [dashed, red curve] and the 
bosonization result in Eq. (21) [solid, orange curve]. 



0.5 while presented data for T < 0.5 are only estimates. Fig- 
ure 8 depicts results for A = 0.5, characteristic results for the 
XY regime. Several conclusions can be drawn directly from 
the obtained results: (a) For T > 1 the decay rate is effec- 
tively T independent, (b) In the regime T < 1 but T > T k 
where FTLM is reliable, it appears that j(q, uo = 0) oc \q\ for 
all T but the T-dependence reveals the vanishing decay rate 
for T — > as expected for low T where the model becomes 
the one of free quasiparticles within the Luttinger ID liquid. 

In the inset of Fig. 8 we show a comparison with the 
bosonization result 7 in Ref. 29 for A = 1 and q = 0, 



7(T) = 7rg 2 T, 



(21) 



where the running coupling constant g is determined by the 
equation 
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We point out that our result, for the smallest possible q\ = 
2tt/L at L = 26, is smaller in the entire T > Tf s regime. It 
should be also noted j(qi,uj = 0) is clearly an upper bound 
to the desired 7(5 — > 0, uj = 0). 

VI. UNIFORM DYNAMICAL SPIN CONDUCTIVITY 

It is by now well accepted that the transport in an integrable 
ID model can be dissipationless at any T > in the "conduct- 
ing" regime. 6,7 For the uniform spin current, given explicitly 
as 

Jq = J ^ (^S r S^_^_ 1 — SrS r +i) , (23) 

r 

one can express the conductivity for q = [cr'(0, uj) = ct'(uj)] 
at T > as 

a'(ui) = 2irDS(u>) + a' res (uj) , (24) 

where the regular part a' Kg (uj) can be expressed in terms of 
eigenstates |n) and eigenenergies e„, 

7r 1 — e~@ u . ^ 
<4 g M = 7 > p n \(n\J^\m)\ 2 5(e n -e m -uj) , 

lj uj 

(25) 

while the dissipationless component with the Drude weight 
(spin stiffness) D is related to matrix elements between de- 
generated states, 

D = J J H Pn\{n\J z \m)\\ (26) 

where p n = e~° Cn jZ are corresponding Boltzmann factors 
and Z is the partition function. 

While there are several analytical as well as numerical re- 
sults supporting D(T > 0) > in the XY regime s = 0, A < 
1 of the integrable AHM, 4 ' 6 ' 10 ' 11 ' 13 there is still controversy on 
the behavior of the regular part a' (ui) in the same regime. 
The upper part of Fig. 9 shows two possible scenarios for 
<j'(uj) differing essentially in the behavior at u) — > 0. Note that 
both scenarios differ also in the integrated normalized spectra 

/(o;) = -L /"dw'ffV), (27) 

— u; 

also presented in Fig. 9, which are much more reliable (mono- 
tonically increasing function) when numerically dealing with 
finite-system results. Some applications of the MF formalism 
combined with a coupling to conserved quantities and pertur- 
bative scattering of nonconserved quantities (similar to the NZ 
approach in this paper) indicate on cr^ eg (0) 7^ 0. 29 ' 31 We will 
in the following present both an analytical argument as well 
as numerical results leading to the conclusion that this is not 
the case. 

As the basis of the reasoning we follow Ref. 32 for the inte- 
grable ID model, characterized by a macroscopic number of 



2ixD (a) 








Figure 9. (Color online) Sketch of two possible scenarios for (a) 
cr'(w) and (b) the normalized integrated spectrum I(ui). 

local conserved quantities Q n , n = 1, . . . , L. Let us consider 
a perturbed Hamiltonian by a fictitious magnetic flux, 19 which 
modifies the exchange (hopping) term in Eq. (1) by the Peierls 
phase factor J — > Jexp(i0). It can be directly verified that 
the perturbed Hamiltonian is still characterized by the same 
conserved quantities, i.e., [H((f>),Q n (</))] = 0. In particular, 
we employ here only Qs(<f)) = Je(4>) representing within the 
integrable AHM the conserved energy current at q = 0. The 
Taylor expansion in 4> of H(4>) and J E (4>) leads to an exact 
relation, 

[if(0), W {<!>)] + [J o m J E (4>)] - 0, (28) 

where W(4>) = d,J E (cj))/dcj) and Jq(4>) = dH((f>)/d(f>. Eval- 
uating the matrix element of the above equation between 
eigenstates |n) and |m), we find 

(n| J o *(0)|m) = (n\W(<j>)\m) ^ ~ ^ , (29) 

where j E are eigenvalues of J E and e n (</>) — e m (<p) = —uj. It 
seems to be plausible 4 ' 32 that j E do not have the same crossing 
points as e n , such that the denominator in Eq. (29) remains 
finite at u> —¥ 0. Taking this into account as well as Eq. (25) 
we see that (n\J§(<fi)\m) oc uj as uj — s- and finally 

lim <t'(uj) oc uj 2 . (30) 

This scenario is depicted in Fig. 9. This behavior of cr'(ui) at 
low frequency is also clearly visible in the normalized inte- 
grated dynamical conductivity I(ui), see lower part of Fig. 9. 
Such a behavior has been already observed in Ref 4, 10, 15, 
and 33. 

To strengthen our arguments we numerically investigate 
(t'(u>) and the normalized integrated I(uj), obtained with ex- 
act (full) diagonalization where (5-peaks are binned in win- 
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0.2 | , , , , 

s = 1/4 A = 0.5 ( a ) 

s = 1/4 A = 1.0 




^ 0.4 - (b) - 
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Figure 10. (Color online) Low-cj dependence of (a) a' ag (to) and (b) 
normalized I(uj). 



dows 5uj = 0.0001. In Fig. 10 we present characteristic re- 
sults on a' (oj) and I(uj) for anisotropics A = 0.5, 1 and dif- 
ferent magnetizations s ~ 0, 1/4. As clearly visible, the regu- 
lar part cr r ' eg (a; — > 0) starts with zero, consistent with Eq. (30). 
However, due to the small system sizes reachable by ED, it is 
hard to differentiate the predicted uj 2 dependence in Eq. (30) 
from some more general power law of a' Kg (uj) oc uj v with 
rj > 0. 

Let us comment on the above results and the form of the 
memory function M(uj) = M(q = Q,u>), as deduced from 
the relation Eq. (8). If one would insert a general form 

7(cj) =70 + oiu r (31) 

into Eq. (8) and assume a nonsingular behavior of M'(uj), it 
is evident that the dissipationless component and <j' ki ,{uj —> 
0) — > are only consistent for 70 = and r > 2. In particu- 
lar, for reproducing the form (30) we need r = 4. Such a result 
clearly puts strong restrictions to proper analytical approaches 
to the dynamical response of integrable systems such as the 
AHM in the "conducting" regime. Still some caution is wel- 
come in the interpretation of our result since both analytical 
argument as well as numerical results should be firm against 
finite-size scaling. In particular, the relation (29) requires the 
uncorrected eigenvalues of e n and beyond the finite size 
uj > 1/L (some our preliminary results confirm this conjec- 
ture) while certain numerical results also reveal low-frequency 
finite-size anomalies 10 although not for presented examples. 



VII. CONCLUSIONS 

Our results confirmed the essential difference of the finite- 
T hydrodynamic behavior between integrable and generic 
nonintegrable quantum many-body systems. For the exam- 



ple of the dynamical spin susceptibility x(Q: w ) an d s P m re- 
laxation function $"(g,w) in the ID generalized AHM we 
showed that a term breaking integrability, in our case non-zero 
next-neighbor repulsion A2 > 0, induces a "normal" diffusive 
behavior in the hydrodynamics, i.e. small (q, uj) regime. This 
is well reflected in the memory-function analysis where the 
spin-current decay rate 7(9, uj) as the central quantity is ef- 
fectively constant, i.e. j(q, uj) ~ 70 in a range of small (q, uj). 
Results are in addition well captured within the perturbation- 
theory approach starting from the XX noninteracting-fermion 
model. 

The integrable model with A2 = we investigated in the 
"conducting" (gapless) regime, i.e. for s = 0, A < 1 and 
s 7^ 0, in this way avoiding the anomalous diffusion as es- 
tablished in the "Mott-insulating" state at s = 0, A > 1. 
From the existence of dissipationless spin transport and finite 
stiffness D(T > 0) > it is evident that the current de- 
cay rate should vanish in the limit (q,oj) — > 0, i.e. j(q — > 
0,uj —> 0) = 0. A nontrivial result of our analysis is, how- 
ever, that in this respect the universality of the noninteracting 
XX model is followed throughout the whole gapless regime 
revealing j(q,oj = 0) cx \q\ at T > 0. While the similarity to 
the XX model occurs in the Iow-cj regime, this is not the case 
at larger u> > u> c = ci- In the latter regime 7(9, w > uj c ) = 
for the XX model while A > induces 7(5, uj) ~ 71 > 
but only weakly (q, a; ) -dependent. In fact, results in this 
uj > uj c regime again resemble the perturbation-treatment 
analysis, which at least within the lowest order clearly fails 
to capture the low-frequency dynamics. 

Nontrivial are also conclusions for the uniform case q = 
where the results for 7(0;) = j(q = 0, uj) are extracted from 
the uniform conductivity cr'(oj). The only consistent possi- 
bility with the dissipationless 6(uj) component and vanishing 
regular part <j'(uj — > 0) = is 7(0; — > 0) oc uj r with r > 2, 
e.g. for r = 4 the analytical argument is correct. This again 
indicates that in an integrable system the current scattering is 
ineffective at low uj at any T in spite of even strong fermion 
interactions at A > 0. 

Although we analyzed only spin-density hydrodynamics, 
one can easily speculate on the behavior of other transport 
properties, in particular the energy density in the context of 
heat diffusion. Since J E is a conserved quantity, the q = dy- 
namical thermal conductivity k(uj) is particularly simple with 
only a dissipationless part and the corresponding thermal- 
current memory function N(uj) = at any T > 0. The 
extension to finite-g energy-density response XE(q, u>) is not 
straightforward, but from the analogy to spin hydrodynam- 
ics and the noninteracting case one can firmly predict that the 
decay rate should vanish again as N"(q,uj = 0) oc \q\. In a 
similar way one can possibly speculate on the T > hydrody- 
namics of other ID integrable models such as the ID Hubbard 
model. 

Experimentally, the most relevant model is the isotropic 
A = 1 ID Heisenberg model realized in several novel 
materials 8 . But in the absence of external magnetic field s = 
case it is also the most controversial one. Our results in Fig. 6 
both for high T ^> 1 as well as for lower T = 0.7 are not 
conclusive due to the restricted system sizes. Nevertheless, 



9 



together with the assumption -f(q — >• 0,0) = (anomalous 
diffusion) one could expect the dependence 7(5, 0) oc q a with 
< a < 1. Since the effective (momentum-dependent) dif- 
fusion coefficient T> oc l/^(q,u> = 0), one could conclude 
that the scaling with the system size L = 2ir/q should follow 
T> oc L a . Indeed at high temperatures such a scaling and di- 
vergence of T> has been recently found in nonequilibriumbath 



scenarios with a ~ 0.5. 
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